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We present a overview of strong correlations in single- wall carbon nanotubes, and an 
introduction to the techniques used to study them theoretically. We concentrate on 
zig-zag nanotubes, although universality dictates that much of the theory can also be 
applied to armchair or chiral nanotubes. We show how interaction effects lead to exotic 
low energy properties and discuss future directions for studies on correlation effects in 
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1. Introduction 

The world of one spatial dimension has long excited theorists studying many body 
physics. The relative importance of interactions as compared to other dimensions 
means that perturbative approaches often fail, and such models often have strongly 
correlated ground states not simply related to the single-particle picture of non- 
interacting electrons. This is coupled with a wealth of analytic techniques, such as 
integrability or Bosonisation which generate non-perturbative results for such one 
dimensional systems, and allow some of the more peculiar properties of strongly 
correlated systems to be investigated. 

Since their discoverjffl in 1991, carbon nanotubes have become an important 
component of modern physics. Aside from their obvious possible applications in 
nanotechnology, they are one of the first example of real one dimensional systems 
where some of the exciting predictions of strong correlations in one dimension can be 
seen experimentally. Indeed, single wall carbon nanotubes (SWNT) are almost ideal 
clean one dimensional systems, with interaction effects dominating over disorder. 

The purpose of this paper is both to give an overview of our current understand- 
ing of strong correlation effects in single wall carbon nanotubes and to provide a 
tutorial introduction to the techniques used to study this for the student who may 
wish to begin work in this area. Obviously, the choice of topics included is highly 
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selective in nature and focuses mainly on recent breakthroughs in our understand- 
ing, for a much more general and complete introduction to carbon nanotubes, we 
refer the reader to one of the many existing books or reviews 033 

In the past decade or so, there have been many papers dealing with 
correlation effects in the ground state of models of single-wall carbon 
nanotubes P I ^ I ^ I ^ I ^ I -'^^ l ^^ l ^^ l ^^ l ^^ l It was realized early on that the low energy physics 
of carbon nanotubes is described by ladder models which are well studied in the 
literature. However, an isolated nanotube has a Coulomb interaction which is only 
very weakly screened - thus adding an extra ingredient to existing theory. Much 
of the early work was based on renomalization group (RG) analysis, but it was 
Nersesyan and Tsvelili!^ in 2003 who first elucidated the underlying strong cou- 
pling phases at half filling in armchair nanotubes by means of a decoupling based 
on the adiabatic approximation and refermionization of the remaining sectors. This 
work was later extended to zigzag nanotube^^^ where the presence of a staggered 
interaction term competes with the umklapp processes to allow the possibility of 
a non-trivial quantum critical point. In the present paper, we review this theory, 
building it up gradually from all the individual processes present in the Hamiltonian 
until we put them together to get a full picture of the ground state of single wall 
carbon nanotubes. 

The structure of the paper is as follows: in Section 2 we look at the single- 
particle picture of electrons hopping on the lattice of a carbon nanotube. In section 

3 we then discuss the important interactions giving rise to strong non-Fermi liquid 
effects. In section 4 we then briefly look at some of the experimental evidence for 
such a theory. We finish with some conclusions and some appendices detailing some 
of the technical and not so technical points involved in the paper. 

2. The non-interacting picture 

Single wall carbon nanotubes are formed by rolling a single layer of graphene into a 
cylinder, along the direction of the superlattice vector (n,m). This wrapping vector 
can be arbitrary, but for the purposes of this review we concentrate on two directions 
with high symmetry - (n,n) armc/iair nanotubes, and (n, —n) zigzag nanotubes. In 
fact, many of the results are universal^^ and so can be applied to other chiral 
nanotubes as well. 

Before looking at nanotubes, one should first understand the band-structure of 
the underlying two dimensional hexagonal graphene lattice.^i^ We start with the 
tight binding model for the two sublattices, a and b 

n = ^tY,ciC^)cbG) + H.c., (1) 
ill) 

where the sum is between nearest neighbors on the hexagonal lattice, 

cj^(?)[cj^(?)] is the creation operator for an electron on sublattice a[b] at lattice site 



November 7, 2008 14:27 WSPC/INSTRUCTION FILE reviewl 



Correlation effects in single-wall carbon nanotubes 3 




Fig. 1. The lattice structure of tlie graphene sheet and armchair (left) and zigzag (right) nan- 
otubes. The underlying lattice is triangular, with a basis of two carbon atoms per unit cell (open 
and closed circles). Note that the two carbon atoms lie at different symmetry positions, but are 
otherwise identical. An (n, m) nanotube is formed by wrapping the graphene sheet along the su- 
perlattice vector ndi + md2 ■ There are two special cases with high symmetry: armchair nanotubes 
when n = m, which in the picture are formed by wrapping the bottom of the sheet to the top; 
and the zigzag nanotubes when n = — m, which are formed by wrapping the sheet left to right. 



i. A Fourier transform gives 

2cos(\/3fc^/2)e-*'=«/2 + e**=«\ f Ca{k)\ 

2cos(V3fc^/2)e*'=«/2 + e-*'^-« ) ycb{k) J 

which on diagonahzation leads to the spectrum 

e(A;^, ky) = ±\/l + 4:COs'^{V3kro/2) + 4:COs{V3k^/2) cos{ky/2). (3) 

This is plotted in Fig. |2] It has the interesting property of being zero only at 
certain points, known as the K points. This means that undoped graphene which 
has precisely one electron per carbon atom so is a half-filled band with chemical 
potential = is a semi-metal, with a zero-dimensional Fermi- surface, also often 
referred to as a zero bandgap semiconductor. This property is due to the symmetry 
between the a and b sublattices and remains, even when we go beyond the tight- 
binding approximation. 

When the graphene sheet is wrapped into a nanotube, the allowed momenta in 
the transverse direction become quantized. The nanotube can then be an insulator 
or metal depending on whether the allowed momenta include the K point or not. 

Much of the early work on nanotubes involved approximating the low-energy 
part of the graphene spectra by Dirac-concs, then quantizing the transverse mo- 
menta to get the allowed energy bands. This approach is very useful for showing 
the universality of conducting nanotubes jl^^l however in this review we will take a 
different approach introduced by LiUjI^l which shows very transparently the relation- 
ship between metallic carbon nanotubes and ladder models. 
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Fig. 2. Tight binding bandstructurc of a single layer of graphene. At half filling, the Fermi surface 
is the discrete set of K points, of which only 2 are independent. 



2.1. Mapping to ladder models 

For the armchair nanotube (n,n), one can make a partial Fourier transform in the 
transverse direction y, 

my 

The tight binding Hamiltonian Eq. [l] then becomes 

m^.q 

We know that the only case when we get gapless modes is when the transverse 
momentum q — 0, which gives the geometry of the ladder shown in FigjsJ^a). 

For the zigzag nanotube, we make an analogous calculation as for the armchair 
nanotube. We make a Fourier transform in the transverse direction, which in this 
case is x. In this case, there is no direct hopping between different states with the 
same ruy, which gives 

'H= [(^\i^y)cq{my + 1) + H.c] (6) 



q my 



where 



with 



t,(m,) = io(g) + (-l)""A(q) (7) 



to{q) = \[t\\+2t^cos{nq/n)\ (8) 
A{q) = i - 2i_L cos(7rg/7i)] 

and q takes integer values j^l < n/2. In the special case where n is divisible by 3, 
the two bands q = ±Q = ±n/3, have a gap A given only by the difference between 
t|| and t±, i.e. curvature effects, which are small and will be discussed later. The 
low energy picture is then one of two chains, illustrated in Fig|3|b). 
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Fig. 3. Top - The underlying ladder models for the low-energy properties of (a) armchair and (b) 
zigzag nanotubes. Bottom - The low energy spectrum of the underlying ladder models. 



2.2. Low energy spectrum and chiral decomposition 

In a one-dimensional system, the Fermi surface is not a surface but a discrete set 
of points. All low-energy physics occurs in the vicinity of these Fermi points. The 
case of the zig-zag nanotube is easiest as the low energy theory consists of two 
decoupled chains. Each of these chains has the usual tight binding spectrum, and 
in the vicinity of the Fermi points this can be linearized so that the spectrum looks 
like two straight lines; one with particles moving to the right at the Fermi velocity 
vf and the other with particles moving to the left at this velocity. 

It is convenient to expand the original lattice operators in terms of smooth fields 
living at one of these two Fermi points - a technique known as chiral decomposition: 

(9) 



where /i is the band-index indicating which 'leg' of the ladder we are sitting on, 
(T is the spin index, and b is the effective lattice spacing. Substituting this into 
the tight-binding Hamiltonian for the ladder model gives the effective Hamiltonian 
Hq = J dxHo, where 



Ho 



-IVp 



(10) 



For the armchair nanotubes, the fact that the ladder is two strongly coupled 
chains means that before we perform the chiral decomposition, we must first diago- 
nalise the Hamiltonian. This is easily done by taking the bonding and antibonding 
combinations of the two legs of the ladder c±^rT = ^ (ci.o- ± C2^a) ■ The spectrum 
is then shown in Fig. [sja) , and clearly shows its relationship to the Dirac cones of 
graphene. 



3. Interactions 

It is well known that reduced dimensions enhance the effects of interactions; in one 
spatial dimension Fermi-liquid breaks down completely, even for arbitrarily small 
interactions and a new paradigm is needed: the Luttinger liquid. We will see how 
such considerations apply to nanotubes. 

It turns out that for undoped isolated nanotubes, there are three important 
interactions to take into account: 
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(1) Unscreened long-range Coulomb interaction, 

(2) Umklapp interactions, 

(3) Dimerization of interactions due to the way we break the C3 symmetry of the 
underlying graphene sheet. 

We will examine each of these interactions in detail in the following sections, to 
identify both their origins in the underlying theory, and the effect they have on the 
low energy electronic properties of single wall carbon nanotubes. The main theo- 
retical tool used to handle interactions in one-dimensional systems is bosonization 
(see Appendix [T|) . However at all stages we supplement the formal arguments with 
cartoons (see Appendix |b]) to aid the reader less familiar with low-dimensional field 
theory. 

As demonstrated in Appendix [T| bosonization is based on the correspondence 

Rna{x) = exp {"iy/n [^^^(a:) - Qf^aix)]} , (U) 

Lf_,a{x) = ^ — exp{~^^/7r[$po■(a:) + Qfj^aix)]} , 

where 3'po-(a;) and Qfj^aix) are a pair of mutually self dual scaler fields with the com- 
mutation relation [^fj^a-{x),Q^a-{y)] = i6^^/6^a-' when y > x and zero otherwise, and 
a is a short-distance cutoff. Basically, this describes the excitation spectrum of the 
one-dimensional Fermi gas as propagating density fluctuations. These fluctuations 
can be categorized by using the linear combinations 

^(*+T+*+i±*-T±*-i). = ^(*+T ±*-T T^"-!), (12) 

which describe total (+) and relative (— ) charge (c) and spin (s) fluctuations. 
The kinetic part of the Hamiltonian then becomes 

nk^n^Y S [{d.eaf + {d.<^af] , (13) 



■±..9±-l 



which is an independent Gaussian model for each of the modes. We will see later 
that this decomposition of modes is very useful in the charge sector, however it is 
somewhat misleading in the spin sector, where the sensible mode decomposition is 
into a singlet and triplet channel, and not the two independent spin-half opjects. 
This is an artifact of using the standard bosonization which is Abelian on a model 
with a non- Abelian SU (2) spin symmetry - however we will see how to fix this later 
and make the symmetries of the model more apparent. 



3.1. Long range Coulomb 

In most situations, we are used to dealing with short range interactions as the 
long-range tail of the Coulomb interaction is rapidly screened. In a nanotube, some 
level of screening can indeed take place through a substrate, however for an isolated 
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nanotube this is weak and one must consider the full long-range nature of the 
Coulomb interaction. 

We begin by looking at what the long-range Coulomb interaction does to the low- 
energy excitations in carbon nanotubes. In terms of the bosonic fields, the density 
p{x) at a point x along the nanotube is given by 



1 



pix) = 



(14) 



and depends only on the total charge sector of the theory, with there being no differ- 
ence between armchair and zigzag nanotubes. The Coulomb interaction is therefore 



H, 



2e' 



Coul 



dx / dy 



d^^+{x)dy^+{y) 



y\ 



(15) 



Now, the Coulomb interaction is eventually screened we make the reasonable as- 
sumption that the screening length Rg is long compared with the radius of nanotube 
R, but short compared with it's overall length L, i.e. i? ^ i?s <C i. Then, assuming 
that the electrons do not form a Wigner crystal, the underlying density is uniform 
along the nanotube and so the above expression can be approximated bj'^ 



Hcoui « ^HRJR) J dx 



(16) 



The sum of the kinetic and Coulomb parts of the Hamiltonian in the total charge 
sector then remains a Gaussian model, but with the fields rescaled - i.e. it is a 
Tomonaga-Luttinger liquid 



dx 



1 



with velocity 




(17) 



(18) 



and Luttingcr liquid parameter — vp/v'^. The long-range nature of the 

Coulomb interaction therefore gives a strong renormalization of the field, giving 
a strong deviation from the non-interacting scaling behaviour of quantities such as 
tunneling density of states. A cartoon picture of this state is given in Appendix |B.1[ 
The Luttinger liquid parameter can also be related to mesoscopic properties of 
the nanotube which experimentally has a finite length.'^ Firstly, there is the level 
spacing 



AE - vp/L, 

and secondly, there is the (capacitive) charging energy 

R 

L R 



(19) 



(20) 
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This means that the Luttinger liquid parameter may be expressed as 

K^{\ + AEj^Ey^'\ (21) 

which allows an independent experimental determination of K to compare with the 
theory. 

Such experiments have been efFectuatecF^ by measuring the tunneling conduc- 
tance. In a Luttinger liquid the tunneling density of states is given by (see Ap- 
pendix A. 2) 

v{uj) ^ \u - ep^ (22) 

where the exponent a depends on the Luttinger liquid parameter, and also whether 
the tunneling is into the end of the middle of the tube: 

abuik^{K-^+K-2)/8 (23) 

aend = (K-^ - l)/4. 

This power law behaviour for the density of states gives similar power law behaviours 
for the transport properties - the conductance at small bias voltage V is G{T) cx T" 
and at large bias voltage we have dl/dV cx F". We will remark more fully on the 
results of the experiments in Section [4] 

As a final remark about the effects of the unscreened Coulomb interaction, the 
strong renormalization in the total charge sector also means that the velocity of 
the plasmons v'^ is much greater than the bare Fermi velocity vp of the excitations 
in the other sectors of the theory. This can be understood by considering that the 
strong interaction gives a very high 'charge stiffness', i.e. changing the total density 
anywhere in the nanotube is difficult because of the long range Coulomb interaction, 
and therefore excitations travel very quickly. This gives rise to a strong separation 
of energy scales between the total charge sector and all other sectors of the theory, 
allowing one to make the adiabatic approximation which will become essential in 
the analysis of the other interactions. 



3.2. Long range Coulomb + umklapp 

The previous section outlined the theory for metallic nanotubes with an unscreened 
density-density Coulomb interaction present. It turns out that such a theory is 
enough for sufficiently doped nanotubes. However an undoped nanotube is exactly 
half-filled and this means one must also consider the physics of Mott insulators - 
in other words one must think about the more detailed structure of the interaction 
arising because the underlying model is actually on a lattice. At half filling, the 
Fermi wavevector kp = 7r/2, which means that terms involving -ikp = 2tt are 
once again uniform (in fact, they are the traditional umklapp terms), and must be 
included. For a cartoon representation of this physics, see appendices |B.2| [R3| and 

El 

We will work with zigzag nanotubes, although the theory as it stands in this 
section is actually more relevant for the armchair case. The reason for this is twofold: 
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firstly it is easier to work with the zigzag nanotubes as the starting point is the ladder 
model with no direct tunneling between the two chains, in contrast to the armchair 
case where the tunneling is stronger than all of the interactions. In fact, it can be 
showrP^ that these two limits are identical within a duality transformation so the 
theory presented here can be equally well applied to the armchair case. So while the 
physics of these two theories is identical, the operator correspondence to the original 
wrapped graphene sheet is not - so we present the results in the form that will be 
useful in later sections when we dig deeper into the physics of zigzag nanotubes. For 
an analogous calculation to this section performed directly on armchair nanotubes, 
see Nersesyan and Tsvelik.'iSl 

One can construct this theory in a number of ways, for example, by applying the 
full structure of the Coulomb interaction on a lattice and proceeding from there, 
or by constructing all allowed terms in the effective Hamiltonian. Here, we will 
use a variation of the latter methodl^iill by constructing the minimal model on the 
original hexagonal lattice which does not contain any spurious symmetries, and then 
using the mapping developed above, compute the low energy effective Hamiltonian. 

It turns out that the minimal model contains both an onsite interaction, U, and 
a nearest neighbour interaction, V. By considering these interactions, along with 
the long range tail of the unscreened Coulomb interaction discussed previously, we 
can construct a complete theory of the low energy physics of armchair single wall 
carbon nanotubes. 





Fig. 4. The umklapp interaction terms. Solid lines represent right moving particles, dashed lines 
left moving particles, fj, is the chain index, and cr the spin index. 



Mapping the U — V model on the original graphene sheet onto the ladder with 
two chains ^ — ±, the interaction is 

ifi ifi 

^^ = ^E^(*H* + 1) (24) 

i 

where n^o- — cjjg-c^cr and n — '^^^n^^r- By carrying out first the chiral decom- 
position on these terms, one obtains a jungle of terms. However, it turns out that 
because of the long range Coulomb interaction, which makes very small, the 
most important terms are the umklapp terms - represented pictorially in Fig. [4] 
This leads to a bosonized Hamiltonian: 



H — Hkin + 'Hcoul + Ti-umk 



(25) 
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Humk = - 2{laY {^^t) 9c COS (y4^<i> 

(33 - 51) COS (\/4^$+) + 53 COS (\/4^e7) - 51 COS (\/4^$; 



The scaling dimension of each of these terms is 1 + , and as is rather small 
from the unscreened Coulomb interaction, they are strongly relevant. The remaining 
terms not written here which are a combination of various backscattering terms all 
have scaling dimension « 2, so are marginal, and will not have any effect on the 
physics. 

Within the U — V model, gc = gi = {U — 3V)/n and g^ = U/n, however we 
will leave gc and gi separate for completeness. The effect of this term was first in- 
vestigated via a renormalization group analysis fHIll however we will take a different 
approactP^ which gives far more insight into the nature of the possible strong cou- 
pling ground states of this model. This approach which centers around an adiabatic 
approximation is based on the strong renormalization of the <i>^ field, and the fact 
that the excitations of this field travel with a velocity much greater than the bare 
Fermi velocity ^ vp- This means that from the point of view of the 'fast' field 
all of the other fields look static, so can be replaced by their instantaneous 
expectation values. The effective problem in the total charge field is therefore a 
sine-Gordon model (see Appendix |A. 3 1 

= Y t^"^)' + + 2na ^^^^ 

where 

U = 5c cos (VATr^-'j+{g3-gi)cos (\/47r$+) -l-gg cos (\/47fe7) -51 cos (\/4^«'7) (27) 

and (3 = The small-/3 sine-Gordon model has a large gap to solitonic 

excitations of the field, but the lowest energy excitations are the breather modes 
Translated into physical terms they are the plasmon's of the original Luttinger 
liquid, but are now gapped. We will come back to the structure of the interaction 
( 26 ) to study the full excitation spectrum of the carbon nanotubes but for now 
we will concentrate on the collective excitations with energies below the gapped 
plasmons. In this region, we can integrate out the <&+ field, which in practice means 
we can replace cos \J 47r<i>+ by its expectation value - note that this approximation 
is also backed up by the RG analysis.'^ This decouples the remaining modes, which 
is clearest to see after refermionization (see Appendix A. 4 1 

3 3 
+ ^ E(-^H^-^' + ^lQ^^l) - J2 XrXI - ^m^xlxl (28) 

j=0 a=l 

where 

T^f = + = + ' 

naVc naVc n 
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■Kavt 



iravtn ' 



(29) 



■navc TraVc n 

The Dirac fermions R f and L f come from the refermionization of the field <i>~ which 
governs vortex excitations. The refermionization of the spin sector in terms of four 
real (Majorana) fermions Xr L' j ~ 0---3 clearly shows the original SU{2) spin 
symmetry is still present, and the excitations separate into a triplet (j — 1 ... 3) 
and a singlet (j = 0) excitation. 



The operators a and /x are Ising operators (see Appendix A.4|, and we use the 
correspondence: 



COs(\/7r(/)+) = criCT2, 
sin(x/7r0+) = /^i^a, 
cos(^/7^07) = CToCTa, 
sin(V7r(/'7) = MoMs, 



COs(v^6'+) ^1(72 

sm{y^6f) = 0-1^2 
cos(^/7^07) = floors 
sin(\/7r6'7) = ctoMs 



(30) 



In general, one has that (ai) 7^ if the corresponding Majorana mass rrii < 0, and 
(fii) 7^ if m, > 0. 

There are eight possible phases depending on the relative signs of m/, mt and 



A 


m/ > 0, ms 


< 0,TOt < 


(cosi 




>^o 




y^O 


(^0) 


y^O 


B 


m/ > 0, TOs 


> 0,TOt < 


(cosi 






(Ma) 


^0 


(^0) 


y^O 


C 


TO/ < 0, TOs 


> 0,TOt < 


(sini 




>^o 


(Ma) 


^0 




y^O 


D 


TO/ < 0, TOs 


> 0,TOt > 


(sini 






(Ma) 


7^0 


(mo) 


^0 


E 


TO/ < 0, TOs 


< 0,TOt > 


(sinI 




>^o 


(^a) 


y^O 


(mo) 




F 


TO/ > 0, TOs 


< 0,TOt > 


(cosi 




>^o 




y^O 


(mo) 




G 


TO/ > 0, TOs 


< 0,TOt > 


(cosi 




>^o 


(^a) 


y^O 


(mo) 


^0 


H 


TO/ < 0, TOs 


> 0,TOt < 


(sini 




>^o 


(Ma) 


^0 


(^0) 


^0 



Possibilities G and H are prohibited if = 5i as would be the case for the 
U-V model of the interactions in nanotubes. These possibilities were all categorized 
by Nersesyan and Tsvelilf^ thus giving the phase diagram (in parameter space) of 
armchair single wall carbon nanotubes (although beware of the difference in sign 
conventions between our derivation for zigzag nanotubes, and their derivation for 
the armchair case). For an unadulterated nanotube one would expect U ^ V > 0, 
i.e. the onsite repulsion is much larger than the off-site repulsion (and both are 
repulsive) , which forces nanotubes to be in our phase F - so we will concentrate on 
only this phase. 

Due to the opposite signs in the masses of the singlet and triplet spin excitation 
modes, there are no local operators on the nanotube that have a non-zero expec- 
tation value in this phase, although there are obviously a choice of non-local order 
parameters. This satisfies the definition of a Haldane spin liquid phase i.e. a gapped 
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phase but with no long range order. The phase can be understood further by defin- 
ing the staggered magnetization as a suitably averaged difference between the local 
spin densities on the two sublattices of the original graphene sheet 

n-(f) = 5„(r)-5fc(r). (31) 

Projecting this onto the low-energy sector of the model gives 

n" - COs(\/7f$7)/io cri/^2Cr3 ^ {cOs{y/Tr^~)) (flo) (ffa)^ '"M ' ^^^"^ 
\criCr2M3/ \A*3/ 

The fact that this contains only one of the Ising models means that the two-point 
correlation function of n~ will display a coherent magnon peak, with a mass gap of 

|TOt|- 

We finally note that as well as the gapped plasmons and the gapped collective 
modes in other channels, the model itself has quasi-particle excitations, which can be 
seen by going back to the original Hamiltonian Eq. [26] A half-soliton simultaneously 
in each of the modes leaves the Hamiltonian invariant, and these are excitations 
carrying charge e and spin 1/2, i.e. they are quasiparticle excitations. We simply 
note that these are gapped on a large energy scale - the most important excitations 
at low energies are the collective spin excitations. 



3.3. Competition of umklapp and single particle gap 

While the previous section should give the full story for armchair nanotubes, which 
have no single particle gap, zigzag nanotubes have a curvature induced gap. As seen 
m Eq. |9] this comes about from a slight difference between iy which goes directly 
along the nanotube and t± which has to partially curve around the waist. It adds 
a term in the Hamiltonian 

i/dhn = A^^(-l)"[ct,(n)c,,^(n+ 1) +i/.c]. (33) 

Now, it was shown by Kane and Mel^^that this gap is small (of order 1/n^) while 
the interaction gaps are all of order 1/n. However, there is another process^"* that 
generates an identical term in the Hamiltonian, but is of order 1/n. This term is 
the staggered interaction, shown schematically in Fig.|5] In fact, being more careful 
about mapping the U — V model of the zigzag nanotube onto the low energy sector, 
one sees that not only the terms in Eq. [24] are generated, but also the staggered 
term 

This staggered term is from the nearest neighbour interaction V , and is entirely 
due to the way the C3 symmetry of the original hexagonal lattice is broken by 
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Fig. 5. The origin of the staggered interaction. When the graphene sheet is wrapped into a zigzag 
nanotube, the interactions do not have a reflection symmetry about any of the lattice points. 

wrapping it into a nanotube - turning the figure on its side shows that no such term 
is present in the armchair case. 

Again, going through the chiral decomposition and bosonization of both the 
exphcit dimerization, Eq. [33] and the dimerized interaction, Eq. |34| it turns out 
they both give the same term in the low energy theorji^l 

'Hdim — g/^O dim 

Cdim = (4/7^a)[cos(^/7^$+)cos(^/7^$;:)cos(^/7f$7)cos(V7^$+) 

+ sin(\/7r<i>+) sin(\/7r«f>;7) sin(\/7r<i>7) sin(\/7r<f>+)]. (35) 

where = V /2n-\-0{l/n'^) is the sum of both the curvature and interaction terms. 
This term must be added to the Kinetic part, the long range Coulomb part and the 
umklapp part to obtain the full Hamiltonian for zigzag nanotubes 

H = TYkin + "^Coul + "^umkl + "^dim (36) 

This Hamiltonian cannot be solved in general, so we consider the physical limit 
U ^ V . The dimerization Tidim may then be considered a perturbation on the 
state created by the umklapp term. As in the previous section, the large U in the 
umklapp term locks the total charge field 0, and one can again perform the 

adiabatic approximation. Noting that (cos(\/7r<f>+)) ^ whereas (sin(-y7r<i>+)) ^ 0, 
one can consider only the first term of the dimerization operator in Eq. |35| which 
after refermionization acquires the following low energy form 

Cdim ~ cos (V7r^'~)/io M 1/^2^*3 ■ (37) 

In the spin liquid phase (see previous section), the operator can be further simplified 
by replacing the operators co&{^/t:<^~) and /ig by their expectation values. The 
resulting operator is: 

= ^cff /^lAi2M3, 

Kq^ P%{cos{^^t)){cos{^<^-)){^lo). (38) 
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All of the physics is now in the spin sector - in some sense, the charge sector of the 
Mott insulating ground state is compatible with the dimerization operator - both 
being happily minimized by the same state. Thus the resulting theory is that of a 
dimerized spin ladder - for a cartoon account of such a system, see Appendix |B.5| 

We see that the dimerization term competes with umklapp processes which 
support the ground state with {(Ji) =/= 0, (/i^) = 0. The cartoon picture shows that 
as 5a is increased to a critical value, an SU{2)i critical point (i.e. the criticality 
of the Heisenberg spin chain) will ensue. Such a claim can be substantiated by 
various theoretical techniques, such as a mapping to a non-linear sigma model-^^or 
a mapping to a generalized Ashkin- Teller model.'^fi' 

The critical point may be estimated as the point where the mass generated by the 
dimerization term alone is equal to the triplet mass. Although both terms originate 
from the nearest neighbour interaction V, they have different scaling dimensions, 
and a simple calculatiorfi^ gives 




where U = Ua/v and V — Ua/v are dimensionless coupling constants. This shows 
that for sufficiently large radius n, the dimerization term will always win, and 
the ground state will be a dimerized phase, a phenomena completely unexpected 
from the scaling of the bare dimerization term. Whether the quantum phase 
transition from the spin liquid to the dimerized phase occurs or not as a function of 
n depends on various non-universal prefactors and is at present an open question. 



3.4. Long range Coulomb + single particle gap 



Another limit of the Hamiltonian Eq. |36] may be solved - when one neglects the 
umklapp processes and considers only the dimerization (either explicit of interac- 
tion induced) and the long range Coulomb terms. It is not clear whether such 
a situation could be realized in nanotubes - but it is nevertheless an instructive 
theoretical exercise in the general properties of the Hamiltonian Eq. 
The total Hamiltonian now is 



36 



n = 



dx 



{d,Q+f + K{d.,^j 



VF_ 

2 



E 

a— {c ,s 



+ (45A/7ra)[cos(\/7i'*i>^) cos(\/7r<i>~) cos(\/7i'<i>7) cos(\/7i'<i>^ ) 

-|-sin(y^$+) sin(V^$^) sin(V^$7) sin(V^$+)] (40) 

One still has the case that the total charge field is much faster than the other 
fields, so the adiabatic approximation may be deployed. However in this case, the 
field is not simply locked at <i>+ = - both the cosines and sines in the dimerization 
term have the same scaling dimension so both must be treated equally. It turns out 
that on integrating out the fast mode, one is left with an 0(6) Gross-Neveau 
modePil, which is integrabld^ and has a spectrum consisting of 3 -I- 6 -t- 3 relativistic 
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massive particles transforming as different representations of 0(6). More details 
about this limit are beyond the scope of this review, we refer the interested reader 
to the original research articlepil 



4. The state of experiments 

Transport experiments on nanotubes have provided some of the best experimental 
data of the Luttinger liquid state to date. The experiments in 1999 of Bockrath 
et a^ clearly show power laws in conductivity as a function of either bias voltage 
or temperature, and all of the data collapses onto a universal scaling curve well 
described by Luttinger liquid theory. Furthermore, as described in Section 3J_ the 
Luttinger liquid parameter can be estimated from mesoscopic effects, and the dif- 
ference in the exponent between tunneling contacts at the end of the nanotube or in 
the bulk of the nanotube can be clearly seen. Although this experiment was actually 
performed on ropes of nanotubes, it was shown beforehancP^ that transport in the 
ropes is dominated by a single metallic nanotube. 

Experiments on isolated single wall carbon nanotubes are more difficult, al- 
though there have been several important studies EMIIMIZI Although these studies 
focussed much more on mesoscopic effects such as Coulomb blockade, they clearly 
show the importance of electron-electron interactions. 

More recently, a spin gap in carbon nanotubes has also been measured for the 
first time by means of an ingenious NMR technique EUlHl In fact, even away from the 
case of half-filling discussed in this paper, interactions are expected to open a spin 
gap.f^ This gap would be much smaller than those predicted by the undoped anal- 
ysis, as they come from marginally relevant backscattering and interband forward 
scattering terms (ignored in this paper) , rather than the strongly relevant umklapp 
ones. We finally comment that to date, neither the Mott insulating phase nor any 
of the unusual properties of the spin liquid phase have been seen experimentally. 



5. Conclusions and future directions 

We have presented a pedagogical overview of the theory of strong correlations in 
metallic single wall carbon nanotubes. We particularly concentrated on the undoped 
case, where the combination of the unscreened Coulomb interaction, umklapp terms 
and a staggered interaction lead to a rich phase diagram - Fig. [6] which is only partly 
understood. 

It is well known that a single chain Mott insulator transforms to a band insu- 
lator when a dimerization term is switched on via not one but two seperate phase 
transitions,!^ one in the spin sector and one in the charge sector. The phase diagram 
of Fig. [6] is reminiscent of this - for small dimerization there is one phase transition 
in the spin-sector at low dimerization, however we know at large dimerization (and 
no Mottness) that the charge is aranged differently. Whether this is reached from 
a crossover or another quantum phase transition is an open theoretical question. 
In fact, there are many open questions in the physics of dimerized ladder models 
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SU(2) QCP 



massive magnoiis 
massless spinons 



massive magnons 



3+6+3 0(6) multiplets 




massive spiiions 



0(6) 



Fig. 



The phase diagram between the known limits 



- different sorts of dimerization may be added (whether relevant for nanotubes or 
not), and it may also be interesting to know how much of a difference screening the 
long range Coulomb interaction would make to this picture. 

Returning to models directly relevant to nanotubes, it is known that the umklapp 
terms (leading to the spin liquid phase) are universal and should be roughly the 
same in all chiral nanotubes as well the as armchair and zigzag cases. However, no 
study has been performed looking at the dimerization of interactions in the general 
chiral case. One might expect chiral to be intermediate between the undimerized 
armchair case, and the maximally dimerized zigzag case, with the chiral angle being 
a plasuible parameter to tune nanotubes to the SU{2)i quantum critical point. 

On the experimental side - nanotubes have already provided some of the nicest 
experimental evidence for the Luttinger liquid state existing in nature. However, 
the predicted Mott insulating phase has never been observed. This may be because 
the nanotubes are naturally doped away from half filling, or possibly the Mott 
insulating state has been seen but the charge gap has been mistaken for one of 
the semiconducting nanotubes. There is as yet no 'smoking gun' experiment to 
distinguish between these two possibilities on an individual nanotube. 

In summary, we have presented an overview of the current (theoretical) under- 
standing of strong correlation effects in carbon nanotubes. With the current rate of 
improvement in experimental manipulation and observation of nanotubes, this is a 
ripe topic with plenty more unanswered questions. 
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Appendix A. Bosonization 

A full account of the bosonization technique is not possible, nor relevant for a review 
of this nature. However, in order to make this article self-contained and to define 
our notations, we present a brief description here. There are a number of good 
reference^^IIMMI for the reader who wishes to delve deeper into this subject. 

The basic idea behind bosonization is that in a model with a linear spectrum, 
excitations consist of particles and holes moving with the same velocity which are 
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therefore coherent. The entire spectrum can thus be reformulated in terms of these 
bosonic particle-hole excitations, i.e. the local density operators, or in other words, 
the plasmons. This is a very useful representation, as the density-density interaction 
remains quadratic, thus allowing the problem to be solved. 
To be more technical, consider the fermionic Hamiltonian 

Hf ^ -ivF{n}d^R - L^d.^L), (A.l) 

which has the single particle propagators (in imaginary time) 

(T,i?(a;,r)i?t(0,0)) = ^^— , {T,L{x,t)LHO,0)) = ^ ^— . (A.2) 
Now, consider the bosonic Gaussian Hamiltonian 

Tis- Y [n2(a;) + (9,ci>(x))2] , (A.3) 
which has the single particle propagator 

{T^x, r)m 0)> - ^ In ( - ^ f ^) + ^ f — ) , (A.4) 

where R is the long distance (low energy) cutoff (i.e. the system size). From this 
expression, we see that at the level of correlation functions, one can decompose the 
field 

^{x,t) = (t)R{x + iT) + (f>L{x - ir), (A. 5) 

i.e. into a 'right moving' part that depends only on the combination x + ir and a 
'left moving' part depending only on a; — ir. A short calculation'^- shows that 

(rY^^il3'i'{x,T)^-tl3'S>{0,0)x ^ g/32(*(a;,T)<E.(0,0)>g-/32(<I.(0,0)<E.(0,0)> ^ / °0 \ ^ M .6) 

\X^ + T^J ' ' 

where uq is a short distance (high energy) cutoff (i.e. the lattice spacing). With this 



expression and the decomposition (A. 5), we see that if we define the correspondence 



i?(x) = ^^e'^'^"^^), L(a;) = ^^e-*^'^^^^), (A.7) 

then the correlation functions between the fermionic model and the bosonic expo- 
nents are identical. In fact, it turn oulP^that the entire spectrum of models A.l and 



[X3| is the same, so the operator mapping |A.7| can be considered an exact equivalence 
between the two models, including when further perturbations such as interactions 
are added to the Hamiltonians. The final ingredient to be added is the correct way 
to calculate the density 

p{x) = R^x)Rix) + L\x)L{x) = -^d^^ix). (A.8) 

It is often useful to introduce the field 



e(x) = -<Pb{x) + <j)L{x) 



(A.9) 



November 7, 2008 14:27 WSPC/INSTRUCTION FILE reviewl 



18 Sam T. Carr 

which satisfies dxQ{x) = Il{x). The fields <i>(x) and Q{x) satisfy the non-local 
commutation relation 



mx),e{y)]^id{x~y). 



(A.IO) 



where 9{x — y) is the Heaviside function. Expressing the Hamiltonian in terms of 
<i>(a;) and Q{x) 



Ti-B 



Vp 



[{dxe{x)f + {dx^{x)f 



(A.ll) 



it is easy to see that Q{x) is the dual field to ^{x), and therefore, the propagator 



(TT-6(a;, t)6(0, 0)) is also given by Eq. A. 4 identical to the $ case 



We now consider perturbing the Hamiltonian |A.1| or |A.3| by a density-density 
interaction 



H' — gp{x)p{x + a). 
In the bosonic case, the Hamiltonian remains quadratic 

I'F 



(A.12) 



Hi 



W = 



V 

2" 



K{d.Mx)f 



1 

K 



(A.13) 



where K = 1/ \/i + g/irvp is known as the Luttinger liquid parameter and v' — 
Vp / K is the renormalized velocity. This equation can therefore be brought back to 
canonical form by the shift of variables <i> = ^/K'^' , = & l\fK. One can then 
calculate fermionic correlation functions in the interacting model 

(T,i?(x,T)i?t(0,0)) ~ (T,e^^[*("'^)-«(^'^)le-^^[*(°^°)-Q(°'°)l) 

_ ^rp gi\/^'I>(a;,T)-'i^7r//fe(a;,T)g-'j\Afi?$(0,0)+i-y7r/A:e(0,0)\j 



x(T^e 



X + IT 
flo 



i(yK+yT7K) 



X + IT \ X^ + T 



X — IT 
(K+l/K~2)/i 



(A.14) 



The Green's function no longer shows a quasi-particle pole at a; = — ir - in a 
Luttinger liquid, all single particle excitations are incoherent. 



A.l. More than one species of Fermion 



The previous section describes the bosonization of a single species of Fermion, how- 
ever in the low-energy effective theory of nanotubes, there are 4 species - 2 spins 
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and 2 bands. The non-interacting Haniiltonian used as a starting point is 

Hf^Y. -^MRlad.R.a - Ll,d,.L^,), (A.15) 

where /i = ± is the band index and cr =t,i is the spin index. Bosonization takes 
place individually for each /icr combination exactly analagous to the previous section 
for the four fields 

R^Ax)^ -^=exp{-iV^[<^>^„x^e^a{x)]} , (A.16) 

where the k^o- are anti-commuting variables (known as Klein factors) which en- 
sure the anti-commutation between the different species of Fermions. In the case 
of nanotubes, these variables have no dynamics, so for most purposes can be safely 
ignored - their only role is to ensure the correct sign of backscattering terms and 
operators representing possible order paramters. We will say no more about them 
in this review, and refer the reader to the original research papers for more details 
about their role in the theory of nanotubes. 

It turns out that instead of the fields <i>^o. which represent density fluctuations 
in the four different fermionic species, the linear combinations 

= ^(ci>+t + a>+x±<i>-T±'i>-i), (A.17) 

which decribe fluctuations in total charge (c+), relative charge (c~), total spin 
(s"*") and relative spin (s^) are a much better basis when interactions are added. 
In this representation, the non-interacting Hamiltonian remains the sum of four 
independent Gaussian models 

nk.n=Y E [{dxOaf + id.<^af] . (A.18) 

a=(c±,s±) 

As an intuitive example, let's see what the single-particle propagator looks like 
when we add the long-range Coulomb interaction - which affects only the total 



charge channel (Eq. 17 in the main text): 



(T,i?+,t(x,r)i?^,^(0,0)) ^n(Tre^^'*"^"'"^"®"^"'"^'e-'\^[*"(°'°)-®»("'"^^ 

a 

where a = ,c'^ , , . The results of the previous section can then be applied to 
each of the four sectors individually as each of the sectors commutes with the other 
ones. The only difference is the coefficient (3 in the exponent e*'^* is now ^7r/4 
rather than ^/tt because of the normalization of the four fields. The result is 

/ X {K+l/K-2)/16 / \ 1/4 / \ 3/4 

(T.i?+,,(x,r)4^(0,0))^(^^) (^^) (^^) '(A-20) 
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the first two terms arising from the c+ sector, and the remaining term being the 
product of the other three sectors of the theory. Similar results arise for each of the 
other possible propagators. 

A. 2. Tunneling Density of States 

The tunneling density of states is calculated as the Fourier transform of the imagi- 
nary part of the retarded Green function 

i^{uj) = — / e"^' 5 G^{x = 0, t)dt. (A.21) 



The retarded Green function is obtained from the Matsubara one (A. 20 1 by the 
analytic continuation t = it + S, so it is seen that 

G^(x = 0, ^ ^-(K+i/K-2)/8+i (A.22) 

is a power law. Now, 

J e''^H-°'dt^ Lo"-^, (A.23) 

so the density of states 

1,{UJ) ^ JK+l/K-2)/S (A 24) 

as stated in the main text. 

So far, everything we have described is for an infinite homogenous system - which 
is a good approximation so long as we stay in the middle of the nanotube. If we 
probe the edge, we must worry about the presence of the boundary. Fortunately 
this is simple within bosonization: no particles can leave, so at the right edge, a 
right mover must be reflected as a left mover, i.e. B){xraax) = L{xmax)- From Eq. 

this implies 



A.9 



A. 201 is modified 



A. 7 we see that this means (t)R{xmax) = (^Lixmax), and from Eq 
that = at the boundary. Therefore, at the boundary Xm, Eq 
to exclude correlations of O: 

(T,i?+.T(x„,r)i?;^^(x™,0)) = [](T.e*\^[*"(-'"^-)le-^%/^*"(-"'°)) 

a 

^ (*t)(i/^-i)/4+i. (A.25) 
which gives the boundary density of states 

y{Lo) ^ (A.26) 



A. 3. The sine- Gordon model 

In the previous sections we showed that density-density interactions in the single 
species case left the bosonic Hamiltonian quadratic, and therefore gave rise to a 
Luttinger liquid, characterized by power law decays of correlation functions. When 
either interband backscattering (e.g. i?^ ^R^^aL^_ ^R+,cr), spin-flip backscattering 
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(e.g. or umklapp terms (e.g. are taken into 

account, this is no longer true, and cosine terms are generated in the Hamiltonian. 

While in the main text, we study the full Hamiltonian for carbon nanotubes 
which include coupling between the channels, here we discuss the basic unit which 
appears in a number of places: the canonical sine-Gordon model in a single channel 

HsG = Y [n'(^) + - cos(/3<i>(x)). (A.27) 

The cosine term would like to lock the field <& at the minimum of the potential 
$(a;) = is g > and <I>(x) = ir/f] if g < 0, however the momentum term lP{x) is 
minimized when ^(x) has large temporal fluctuations. The relevant importance of 
each of these terms is controlled by the parameter /?. The scaling dimension of the 



operator cos (/?<!> (a:)) can be read off from Eq. A. 6 d = /3^/47r. When the scaling 
dimension is greater than 2, correlation functions of this operator do not diverge 
(at low energy) when integrated over position and time. A scaling dimension of less 
than 2 indicates a divergence will occur when perturbation theory is applied - and 
so the operator may substantially change the ground state. 

The sine-Gordon model is integrable, and is one of the most studied models 
in physics.fS^ Here, we will simply give a summary of the main results without 
derivation. If > Stt, then the momentum term wins, and the cosine perturbation 
fails to lock the field ^{x) - in technical terms we say it is an irrelevant operator. In 
this case, we can forget about it's presence and deal only with the Gaussian model. 

If however, < Stt, then the cosine term wins, and the field ^{x) is 
locked, giving a non-zero expectation value to (cos(/3/2<I>(x))) if g > and 
(sin(/3/2<I>(a;))) 7^ if g < 0. There is a spectral gap to excitations, and the exci- 
tations are given by solitons in the field <i>(x), interpolating between two adjacent 
minima of the cosine potential. These solitons can be visualized as in the cartoon 
Appendix. |B.3[ An anti-soliton is then the mirror image of a soliton. 

Furthermore, there is a residual interaction between solitons and anti-solitons - 
which is repulsive if 47r < < Stt and attractive if < in. In the attractive case, 
the solitons and anti-solitons can form bound states (with an energy gap less than 
the individual solitons), which are once more similar to the density waves of the 
original Gaussian model - this is what happens in the total charge channel on the 
half-filled nanotubes (see main text) . At exactly — 4tt, the solitons are free and 
may once more be re-expressed in term of fermions. 

A. 4. Refermionization and Majorana fermions 

Refermionization is due to the equivalence of a model of massive Dirac fermions to 
the = 47r sine-Gordon model, demonstrated by taking the Hamiltonian 

n,G-F = -ivpiR^'d^R - L^d^L) + ig{R^ L - R) (A.28) 

and bosonizing according to the rules ] A. 7] The result is the sine-Gordon Hamiltonian 



A.27 (the cosine is obtained rather than a sine from the correct treatment of the 
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product of the exponents of two non-commuting operators, <i> and Qj^. In this 
case, the bosonization can work backwards and is known as refermionization. The 
Fermions represent the non-interacting soHtons of the sine-Gordon model. It must 
however be remembered that the sine-Gordon model usually appears in one of the 
sectors of the theory invoked by the linear combination |A.18| so are not simply 
related to the original Fermions i?^^^ of the model (in fact the relationship is in 
general not even local). 

When treating the spin sector of nanotubes, it is convenient to re-express the 
Dirac Fermions in terms of Majorana (real) fermions: 

1 _ i?^ + i? 2 _R'^ - R 
Xr- 2 ' 2i ' 

Xl-—, Xl-—. (A.29) 
The Hamiltonian of the Majorana fermions is then the kinetic part: 

2 

Um = -IV F (Xr^^Xr - Xld.xl) , (A.30) 
and the possible mass terms are 

XrxI + Xlxl = ^R^L - L^R) = cos{^<^) 

XrxI - xlxl = R^L^ + LR= cos(V4^e) 

(A.31) 

In the main text, the refermionization of the spin sector in terms of Majorana 
fermions comes about by carrying out this procedure for the two sectors and 
<I>~, thus arriving at 4 Majorana fermions. 

The utility of this representation comes about from the equivalence between a 
single species of Majorana fermion and the transverse field Ising model 

Hi = -jY.<yt<yUi + hJ2<yf (A.32) 

i i 

after carrying out a Jordan- Wigner transformation EIEI! The mass of the Majorana 
comes out to be to = |J| — \h\, with the Z2 quantum critical point at |J| = \h\. If 
I J| > \h\, then this model is in it's ordered phase (cr^) 7^ 0, whereas in the opposite 
case \h\ > \ J\, the so called disorder parameter /if — Ylj^^ cTj is non-zero {fi^) ^ 0. 
By bosonizing two copies of the Ising model, one can relate cr and fi to the bosonic 
fields <i> and bEIIMI For details of this, see the literature - the correspodence we 
use is shown in the main text, and allows the strong coupling phases of carbon 
nanotubes to be understood. 



Appendix B. Cartoon pictures of the physics of Mott insulators, Spin 
chains and Spin ladders 



November 7, 2008 14:27 WSPC/INSTRUCTION FILE reviewl 



Correlation effects in single-wall carbon nanotubes 23 

Here, we provide a strong coupling 'cartoon' picture of the physics of the correlated 
states of carbon nanotubes. These are provided both as a way to understand the 
one-dimenionsal physics involved without delving into the details of bosonization, 
and as a way of visualizing the underlying field theory. We stress that these are 
cartoons and are no substitute for a proper calculation, however we feel that they 
give a very intuitive window into the physics of one dimensional systems, which is 
sometimes lost in the technicalities of field theory. 

B.l. One dimensional interacting electrons 

First, consider electrons confined to move in only one dimension, i,e. on a line. 
Even with arbitrarily small interactions, the electrons cannot go through each other, 
therefore there are no (coherent) single particle electron like excitations, and all 
excitations are plasmons. 

• cT # ^ m % 

This is a basic picture of the Luttinger liquid state. 
B.2. Half filling in a lattice - the Matt insulator 

Now consider a one-dimensional lattice model - the Hubbard model with an energy 
cost U for two particles to be on the same site, and a nearest neighbor hopping 
integral t. At half filling, the electrons want to avoid the energy cost [/, so have 
exactly one electron per site - the Mott insulator. It turns out that in contrast 
to higher dimensions, a one dimensional system will form a Mott insulator at half 
filling for arbitrarily small on-site interaction U . 

Hubbard model: ^ ^ ^ _^ 

u 

Half filling: 

--(^--^--^--^--(^--^--^--(^-- 

There is a gap of order U to any charged excitation, but there is a residual anti- 
ferromagnet exchange interaction of order /U from virtual hopping terms. All of 
the interesting low-energy physics is in the spin sector - so we need to consider the 
physics of one dimensional spin chains. 

B.3. Spin Chain - deconfined solitons 

Consider making an excitation in an antiferromagnetic spin chain by fiipping the 
spin in the center - this changes the spin of the system by 1 and creates two un- 
satisfied bonds - as shown in figure (a). These unsatisfied bonds however can move 
independently, splitting into two solitonic excitations known as spinous, each car- 
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rying a spin 1/2 - figure (b). 

We must stress here again that this is a cartoon picture - in fact, the Heisenberg 
spin chain has only quasi-long-range order (i.e. power law correlations), even at zero 
temperature, and so never actually forms the Neel state pictured. Nevertheless, the 
basic physics is represented by the cartoon is correct. 



B.4. Spin Ladder - confinement of solitons 

At this stage, we must remember that the low-energy model relevant for nanotubes 
is actually a ladder model, and not a single chain. While this makes very little 
difference in the charge sector - at half filling the ladder is still a Mott insulator, 
the physics of the spin ladder turns out to be very different from the spin chain. 



-fitifffffHifif 
-IfHifH-lHfHi- 

Suppose two solitons form on one of the legs of the ladder and start to move apart. 
Then all of the bonds from the first leg to the second between the two solitons 
are unsatisfied - thus as the solitons move further apart, the energy increases lin- 
early creating a confining potential (actually, very similar in style to the confining 
potential between two quarks in quantum chromodynamics) . The excitations can 
then roughly be split into two classes: either the two solitons pictured must move 
together, which is a very rough cartoon of the triplet excitations, or a soliton on 
one chain moves coherently with an antisoliton on the other chain - which is very 
roughly the singlet excitation. It turns out that the confinement energy leads to 
both of these classes of excitations being gapped - although again one must remem- 
ber that in reality, there is only quasi-long-range order and the Neel state is not 
truly formed. 



B.5. Dimerized spin ladder - quantum criticality 

Finally, we consider the dimerized spin ladder important for the zig-zag nanotubes. 
In the absence of dimerization, the model is that of a uniform spin ladder, as pictured 
in section B.4 which has gapped triplet and singlet excitations (or in other words. 



excitations consisting of confined spinons). However, now turn on the dimerization 
and consider tuning it very carefully as shown in the figure. 
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A = J|| 
J I = 2J 



At the critical point, the ladder becomes equivalent to a single chain - |B.3| - which 
has gapless excitations consisting of deconfined solitons. This is the so called SU {2)i 
quantum critical point expected in zig-zag nanotubes. Note that a full calculation 
shows that the critical point is in fact not just a point as in the cartoon, but an 
entire line Ac( Jii, Jl). 



References 

1. S.Ijima, Nature (London) 354, 56 (1991). 

2. R.Saito, G.Dresselhaus and M.S.Dresselhaus, Physics properties of carbon nanotubes 
(Imperial College Press, London, 1998). 

3. J.-C. Charlier, X.Blase and S.Roche, Rev. Mod. Phys. 79, 677 (2007). 

4. C.Dekker, Phys. Today 52(5), 22 (1999). 

5. R.Egger and A.O.Gogolin, Phys. Rev. Lett. 79, 5082 (1997). 

6. C.Kane, L.Balents and M.P.A.Fisher, Phys. Rev. Lett. 79, 5086 (1997). 

7. L.Balents and M.P.A.Fisher, Phys. Rev. B 55, Rll 973 (1997). 

8. Yu.A.Krotov, D.-H.Lee and S.G.Louie, Phys. Rev. Lett 78, 4245 (1997). 

9. H-H. Lin, Phys. Rev. B 58, 4963 (1998) 

10. R.Egger and A.O.Gogolin, Eur Phys Lett. (1998). 

11. H.Yoshioka and A.A.Odintsov, Phys. Rev. Lett. 82, 374 (1999). 

12. A.Odintsov and H.Yoshioka, Phys. Rev. B 59, RIO 457 (1999). 

13. A.A.Nersesyan and A.M.Tsvelik, Phys. Rev. B 68, 235419 (2003). 

14. S.T.Carr, A.A.Gogolin and A.A.Nersesyan, Phys. Rev. B 76, 245121 (2007). 

15. P.R.Wallace, Phys. Rev. 71, 622 (1946) 

16. M.Bockrath, D.Cobden, J.Lu, A.Rinzler, R.Smalley, L.Balents and P.L.McEuen, Na- 
ture 397, 598 (1999). 

17. C.Wu, W.V.Liu and E.Fradkin, Phys. Rev. B 68, 115104 (2003). 

18. C.L.Kane and E.J.Mele, Phys. Rev. Lett. 78, 1932 (1997). 

19. M.A.Martin-Delgado, R.Shankar, and G.Sierra, Phys. Rev. Lett. 77, 3443 (1996). 

20. Y.J. Wang and A.A.Nersesyan, Nucl. Phys. B 583, 671 (2000). 

21. L.S.Levitov and A.M.Tsvehk, Phys. Rev. Lett. 90, 016401 (2003). 

22. A.B.Zamolodchikov and Al.B.Zamolodchikov, Ann. Phys. (N.Y.) 120, 253 (1979). 

23. M.Bockrath et al, Science 275, 1922 (1997). 

24. S.J. Tans et al, Nature (London) 386, 474 (1997). 

25. S.J. Tans et al, Nature (London) 394, 761 (1998). 

26. D.H.Cobden et al, Phys. Rev. Lett. 81, 681 (1998). 

27. W.Liang et al, Phys. Rev. Lett. 88, 126801 (2002). 

28. P.M. Singer, P.Wzietek, H.AUoul, F.simon and H.Kuzmany, Phy. Rev. Lett 95, 236403 
(2005). 

29. B.Dora, M.Gulacsi, F.Simon and H.Kuzmany, con-mat/0707.3348 (2007). 

30. M.Fabrizio, A.O.Gogolin, and A.A.Nersesyan, Phys. Rev. Lett. 83, 2014 (1999). 

31. A.O.Gogolin, A.A.Nersesyan, and A.M. Tsvelik, Bosonization and Strongly Correlated 
Systems, (Cambridge University Press, Cambridge, England, 1998). 

32. T.Giamarchi, Quantum Physics in One Dimension, (Oxford University Press, Oxford, 
England, 2004). 



November 7, 2008 14:27 WSPC/INSTRUCTION FILE reviewl 



26 Sam T. Carr 

33. D. Senechal, |cond-mat/9908262| (1999). 

34. G.E. Japaridze, A.A.Nersesyan and P.B.Wiegmann, Nucl. Phys. B 230, 10 (1984). 

35. D.G.Shelton, A.A.Nersesyan and A.M.Tsvelik, Phys. Rev. B 53, 8521 (1996). 



